•I  ,-\ 
■i ' 

I  vi 

■'jJ  : 
1  - 


i’t'TBfrtfaMj'f-* 


AD-756  832 

EXTREMES  OF  HYDROMETEORS  AT  ALTITUDE 
FOR  MIL-STD-2I0B.  SUPPLEMENT:  DROP- 
SIZE  DISTRIBUTIONS 

Paul  Tattleman,  et  al 

Air  Force  Cambridge  Research  Laboratories 
L.  G.  Hanscom  Field,  Massachusetts 

20  October  1972 


DISTRIBUTED  BY: 


National  Technical  information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va,  22151 


h.  V  •.  '*■' 


.^:,;.‘;??-s- 


MR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 


L.  G.  HAMSCOM  FIELD,  BEDFORD,  MASSACHUSETTS 


MIL-STD-210B 


Supplement  -^Drop-Size  Distributions 


PAUL  TATTELMAN 
NORMAN  SISSENWINE 


Reproduced  by 

NATIONAL  TECHNICAL 


INFORMATION  SERVICE 

U  S  Department  of  Commerce 
Spr iny f If* I H  VA  22151 


ASP  FORCE  SYSTEMS  COMMAND 

United  States  Air  Force 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  RP*D 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 

t.  ORIGINATING  ACTIVITY  (Cotjkfate  author) 

Air  Force  Cambridge  Research  Laboratories  (LK) 

L.  G.  Hanscom  Field 

Bedford,  Massachusetts  01730 

20.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 
to,  group 

J.  REPORT  TITLE 

EXTREMES  OF  HYDROMETEORS  AT  ALTITUDE  FOR  MIL-STD-210B 

Supplement  -Drop-Size  Distributions 

«.  OEICRIPTIVE  NOTE*  (Type  of  report  tnd  inclusive  dales) 

Scientific,  Interim. 

j.  auyhorMI  (first  none,  middle  initial,  tost  none) 

Paul  Tattelman 

Norman  Sissenwine 

ft.  REPORT  DATE 

20  October  1972 

70.  TOTAL  NO.  OF  PAGES  76.  NO.  OF  REFS 

IS  8 

to.  CONTRACT  ON  GRANT  NO. 

6.  PROJECT,  TASK.  WORK  UNIT  NO*.  8624  01  01 

C.  000  ELEMENT  62  10  IF 

d.  POO  SUBELEMENT  681000 

_ 

SO.  ORIGINATOR'S  REPORT  HUMBERTS; 

AFCRL-TR-73-0008 

tb,  otherreporT  HOiS)  (Any  other  numbers  that  may  be 
o .signed  Ibis  report)  APSQ--Nm~21?3 

».  distribution  statement 

Approved  for  public  release;  distribution  unlimited* 

It.  SUPPLEMENTARY  NOTES 

TECH,  OTHER 

ta.  sponsoring  military  activity 

Air  Force  Cambridge  Research 
Laboratories  (LK) 

L.  G.  Hanscom  Field 

Bedford,  Massachusetts  01730 

IS.  ABSTRACT 


In  the  basic  report  for  this  supplement,  estimates  were  provided  for 
extreme  rainfall  intensities  and  associated  liquid  water  content  at  altitude 
with  0.  1  and  0.  5  percent  probabilities  in  the  most  severe  month  in  the  rainy 
tropics.  Analagous  values  were  also  included  for  two  all-time  record  rain 
storms:  one  for  1  min,  and  the  other  lasting  42  min.  In  this  brief  supplement, 
drop-size  distributions  are  developed  for  the  1-min  record,  and  the  0.  1  and 
0.  5  percent  "Worst  Month  Tropics"  from  a  theoretical  expression  that  best 
fits  the  available  empirical  size  distributions.  Liquid  water  content  obtained 
from  an  integral  of  this  theoretical  expression  is  compared  to  the  liquid  water 
content  obtained  by  arithmetically  adding  the  volumes  of  each  drop  for  three 
drop-size  distributions. 
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Extremes  of  Hydrometeors  at  Altitude 

for  MIL-STD-210B 

Supplement  — Drop-Size  Distributions 


1.  INTRODUCTION 

Sissenwinc  (1972),  in  the  basic  report  for  this  supplement,  provided  estimates 
of  extreme  rainfall  intensities  and  associated  liquid  water  content  at  altitude  for 
0.  5  and  0.  1  percent  probabilities  in  the  most  severe  month  in  the  rainy  tropics. 
Analagous  values  were  also  included  for  the  most  severe  all-time  recorded 
42-nun  and  1-min  rain  storms.  These  values  were  proposed  for  inclusion  in 
MIL, -STD-2 10B,  "Climatic  Extremes  for  Military  Equipment.  "  At  a  subsequent 
DoD  Task  Group  meeting  on  this  DoD  document,  a  need  for  drop  sizes  and  numbers 
was  expressed.  In  this  brief  addendum,  drop -size  distributions  by  altitude  are 
developed  for  the  1-min  record  (highest  ever  recorded)  and  the  0,  1  and  0.  5  per¬ 
cent  "Worst  Month  Tropics"  profiles.  Liquid  water  contents  are  computed  for  a 
few  of  these  drop-size  distributions  by  arithmetically  adding  the  volumes  for 
each  drop.  These  are  compared  to  estimates  of  liquid  water  content  from  an  ex¬ 
pression  provided  in  the  basic  report  which  provides  integrated  water  mass. 


2.  DROP-SIZE  DISTRIBUTION 

Dyer's  (1969)  review  and  summary  of  findings  on  drop-size  distributions 
serves  as  a  basis  for  developing  specific  distributions  herein.  From  an 

(Received  for  publication  20  October  1972) 
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investigation  of  summer  rain  over  Ottawa,  Canada  with  intensities  of  1  to  23  mm/hr, 
Marshall  and  Palmer  (1948)  found  that  except  for  small  diameters,  the  size  distri¬ 
bution  of  raindrops  can  be  estimated  using  an  exponential  expression  of  the  form 

Nd  -  Noe'AD  (1) 

where  is  the  number  of  drops  of  unit  size  range  per  unit  volume,  D  is  the  drop 
diameter,  Nq  is  a  fictitious  value  of  Nq,  the  intercept  with  the  ordinate  when  D  is 
zero  (empirical  studies  indicate  a  cutoff  in  raindrop  size  betweer  0  and  0.  5  mm), 
and  A  is  a  variable  which  is  dependent  upon  the  intensity  of  the  precipitation.  The 
relationship  of  A  with  rain  intensity  is  fixed  by  the  expression 

A  =  aRb  (2) 


where  R  is  the  intensity  and  a  and  b  are  constants.  These  values  vary  with  the  in¬ 
vestigator.  Dyer  also  indicates  that  the  size  of  the  median  volume  raindrops  (that 
is,  the  diameter  that  divides  drops  of  larger  and  smaller  diameter  into  groups  of 
equal  total  vomme),  DQ,  can  be  expressed  in  millimeters  by 


D„ 


3,  67 
A 


(3) 


where  A  (mrr.  ~A)  is  dependent  on  the  precipitation  intensity,  R  (mm/hr). 

Several  authors  have  derived  values  for  Nq  and  A  from  empirical  data  for 
various  locations  and  precipitation  characteristics.  Dyer  (1969)  summarizes  the 
results  of  most  of  these.  Since  concern  herein  Is  with  high  intensity  rainfalls,  it 
would  seem  logical  to  use  a  distribution  derived  from  thunderstorm  data,  but  there 
is  some  diversity  among  investigators  of  drop-size  distributions  on  the  expres¬ 
sion  that  best  fits  such  distributions. 

Mueller  and  Sims  (1970)  have  recently  presented  empirical  spectrums  of 
drop-size  distributions  near  the  surface  which  were  obtained  by  photographic 
observations  for  various  rainfall  intensities  in  the  tropics,  including  higher  rates 
of  fall  of  greatest  importance  in  military  design.  These  are  given  in  Figure  1. 

To  this  figure,  theoretical  drop-size  distributions  have  been  added  for  a  rainfall 
Intensity  of  40  mm/hr.  These  were  computed  from  Eq.  (1)  for  NQ  and  A  values 
suggested  by  three  investigators: 


Source 

Type  of 
Precipitation 

No 

(mm*  *m"2) 

A 

(mm’1) 

Do 

£nm) 

Marshall  and  Palmer  (1948) 

Ottawa,  Summer 

8,  000 

4.1^’ 

0.9R~^2i 

Jones  (1956) 

Thunderstorm 

3BBR1’02 

2.  5R"0,  05 

*—» 

CO 

O 

o 

Joss  et  al  (1970) 

Thunderstorm 

1,400 

3.  OR’0,21 

1.  22R0,  21 

NUMBER  OF  DROPS  /m*  0.1mm 
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DIAMETER,  mm 

Figure*  1.  Avorage  Drop'  si?,o  Spectra  for  Natural  Rainfall  Rates  Occurring  0,  01, 
0, 1,  0.  5,  and  1,  0  Porcont  of  the  Time  (Dashed  Linos),  These  curves  are  for 
tropical  climates  (Muollor  and  Sims,  1070),  The  theoretical  distributions, 
continuous  linos,  for  the  throe  sources  shown  wore  computed  for  a  precipitation 
rate  of  40  mm/hr 
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Computations  were  made  for  each  1-mm  class  interval  starting  with  0,  5  to 
1.4  mm.  An  examination  of  Figure  1  indicates  that  the  distributions  provided  by 
the  first  two  investigators  conform  best  to  Mueller  and  Sims'  data.  The  Joss 
distribution  tends  toward  an  unrealistically  high  number  of  very  large  drops  while 
providing  a  relatively  low  estimate  of  the  number  of  small  drops.  The  straight 
lines  of  the  theoretical  distributions  of  Eq.  (1)  do  not  agree  with  curves  of  the 
empirical  data  below  1.  5  to  2.0  mm  diameter.  Tt  is  apparent,  however,  that 
Eq.  (1)  overestimates  the  number  of  smaller  drops  since  it  indicates  the  greatest 
number  when  the  diameter  is  zero. 

The  median  size  raindrops,  Dq,  for  rainfall  rates  of  48,  188,  and  1860 
mm/hr  are: 

Intensity  (mm/hr) 


48 

188 

1860 

(mm) 

(mm) 

(mm) 

Marshall  and  Palmer  (1948) 

2.  03 

2.  70 

4.  37 

Jones  (1956) 

1.  80 

1.  92 

2.  16 

Joss  et  al  (1970) 

2.  80 

3.  66 

5.  93 

Comparing  these,  the  Jones  (1956)  medians  are  least  sensitive  to  rain  in¬ 
tensity,  remaining  close  to  2  mm.  This  would  indicate  that  a  heavier  rainfall  has 
more  drops  of  all  sizes  which  is  in  agreement  with  the  curves  of  empirical  dis¬ 
tributions  of  Mueller  and  Sims  (1970)  in  Figure  1,  except  for  a  minor  discrepancy 
in  the  two  heaviest  intensities  for  sizes  below  about  1.  1  mm  .  ;i.>  :•  two 

expressions  provide  unrealistically  large  medians  at  the  h;...'.  ■••tt.  r,  ror -re  1 

shows  that  Jones'  values  for  N  and  Aare  most  suitable  «<■:  s.h  utt»me  *  a  •> .fa  11 
rates  with  which  wo  are  concerned,  and  utilizing  them  is  n.u'teni  wait  tin-  basic 
report.  The  resulting  equation  for  the  distribution  of  drop  ‘  i/.ci;  iw 

Nd  *  389R1,  02  exp  (4) 

u  o 

where 

D0  «  1.48R0,05  .  (B) 

Equation  (4)  can  be  Integrated  to  express  the  total  precipitation  water  content 
in  tile  form  of  tho  equation 

M-O.Q52R0, 87  (6) 

where  M  is  the  water  contont  (g/rn  )  and  R  the  intensity  (mm/hr).  Slseenwlne 
(1972)  utilized  this  equation  in  developing  his  precipitation  profiles. 


Sissenwine  (1972)  considers  the  precipitation  in  the  extreme  storms  presented 
in  the  basic  study  to  be  all  liquid  below  4,  5  km.  Between  4.  5  and  10  km  these 
storms  contain  ice  and  liquid  particles  mixed,  becoming  mostly  solid  (snow)  above 
10  km.  Diameters  of  solid  particles  can  be  computed  with  equations  presented  by 
Dyer  (1969)  which  are  somewhat  more  complex  than  for  liquid  precipitation.  The 
hypothetical  aspects  of  applying  these  equations  to  extreme  altitudes  and  intensities 
beyond  that  from  which  they  were  developed,  however,  does  not  warrant  this  addi¬ 
tional  refinement.  Therefore,  equivalent  liquid  water  using  Eqs.  (4)  and  (5)  will 
be  assumed  for  higher  altitudes. 

Precipitation  rates  at  the  surface,  1  km,  2  km,  and  higher  2 -km  intervals  up 
to  18  km  were  interpolated  for  the  proposed  MIJ,-STD-210B  from  Table  3  of 
Sissenwine  (1972).  Median  diameters  and  the  distributions  of  drop  diameters, 
computed  for  these  levels  with  Eqs.  (4)  and  (5),  are  shown  in  Tables  1,  2  and  3. 

To  obtain  the  drop-size  distributions  for  1-mra  intervals  in  Tables  1,  2  and 
3,  the  midvalue  for  each  millimeter  class  interval  of  the  range  of  drop  diameters 
was  used  for  D  (for  example,  1  mm  was  used  lor  D  to  comp  ite  the  number  of  drops 
in  the  0.  50  to  1.  49-mm  class  interval,  etc. ), 

As  already  noted,  Eq.  4  indicates  that  as  the  drop  size  D  approaches  zero, 

1  02 

the  number  of  drops  continues  to  increase  toward  339R  ’  .  This  is  unrealistic, 

since  it  leads  to  the  largest  number  of  drops  at  zero  diameter.  This  is  a  result 
of  fitting  the  distribution  of  raindrop  sizes  with  a  simple  e  .u.  tion  which  neglects 
physical  limitations  that  govern  both  large  and  small  sizes  attainable  in  nature. 

The  most  important  contribution  to  the  total  water  content,  however,  comes  from 
the  center  of  the  distribution  at  diameters  of  2  to  about  5  mm.  The  assumption 
fhat  the  logarithm  of  the  number  of  drops  varies  inversely  with  the  size,  as  in 
Eq.  (4),  is  borne  out  by  Figure  1.  Figure  1,  however,  also  shows  that  this 
theoretical  distribution  departs  from  the  shape  of  the  empirical  curves  at  diameters 
less  than  2  mm  when  Jones'  (1956)  values  for  N  and  A  are  used.  This  short¬ 
coming  should  not  be  of  great  importance  to  the  total  liquid  water  content  since 
the  difference  in  water  mass  for  these  very  small  drops  is  small. 

Another  important  input  from  Figure  1  is  that  the  empirical  data,  in  addition 
to  showing  that  the  number  of  drops  has  a  maximum  at  sizes  of  1  to  1.  5  mm 
depending  upon  intensity,  also  indicates  there  are  no  drops  of  less  than  roughly 
0.  5 -mm  diameter.  A  sharp  cutoff  at  0,  5  mm  is  questionable,  but  the  relatively 
slower  terminal  velocity  of  the  smallest  drops  will  lead  to  collision  and  coalescence 
with  faster  falling  larger  drops,  providing  a  sharp  diminution  of  drops  at  some 
diameter  evidently  not  too  much  lower  than  0.  5  mm.  For  this  reason,  the  smallest 
size  class  used  for  the  distributions  provided  in  Tables  1,  2  and  3  is  0.  5  to  1.  4  mm. 
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Table  1.  Number  of  Drops  (Liquid  or  Equivalent  Liquid)  per  r 
1-min  World-Record  Rainfall 


Altitude 

(km) 

Rainfall 

Median 

Diameter  (mm) 

(mm/min) 

(mm) 

0.  5-1.4 

1.5-2.  4 

2.  5-3.4 

3.  5-4.4 

4.5-5 

0 

31 

2.2 

158.624 

29.915 

5642 

1064 

201 

1 

34 

2.2 

174,297 

32,870 

6199 

1169 

220 

2 

36 

2.2 

184,761 

34,844 

6571 

1238 

234 

4 

41 

2.2 

210,970 

39,787 

7503 

1415 

267 

6 

42 

2.2 

216,220 

40,777 

7690 

1450 

274 

8 

30 

2.2 

153,407 

28,930 

5456 

1026 

194 

10 

20 

2. 1 

93,699 

16,321 

2843 

495 

86 

12 

14 

2.  1 

65.  123 

11,  344 

1976 

344 

60 

14 

9 

2.0 

38,024 

6,069 

669 

155 

25 

16 

4 

1.9 

15,097 

2,  188 

317 

46 

7 

18 

1 

i.e 

3,298 

428 

56 

7 

1 

20 

0 

- 

- 

- 

- 

- 

Table  2.  Number  of  Drops  (Liquid  or  Equivalent  Liquid)  per  m"  Equalled 
or  Exceeded  With  0.  1  Percent  Probability,  Worst  Month,  Wet  Tropics 


Rainfall 

Median 

Diameter  (mm) 

_ 

(km) 

(mm/min) 

(mm) 

0.  5-1.4 

I. 5-2. 4 

2. 5-3. 4 

3.  5-4.4 

4.5-5. 

ntwni 

0 

3.  13 

1.9 

11.755 

1,704 

247 

30 

g, 

l 

1 

3.38 

1.9 

12,714 

1,843 

207 

39 

6 

1 

2 

3.  63 

1.  9 

13, 674 

1,982 

287 

42 

6 

1 

4 

4.  11 

1.9 

15,621 

2,249 

326 

47 

7 

1 

6 

4.23 

1.  9 

15, 983 

2,310 

330 

49 

7 

l 

0 

2.97 

1.9 

11,  143 

1,815 

234 

34 

5 

l 

10 

2.03 

1.9 

7,559 

1,005 

159 

23 

3 

<  1 

12 

1.  38 

l.  8 

4,580 

606 

78 

10 

1 

<  l 

14 

0.  85 

1.  B 

2,794 

364 

47 

0 

1 

<  1 

16 

0.43 

1.7 

1,237 

143 

16 

2 

<  1 

<  l 

18 

0.  00 

1.6 

210 

22 

2 

<  1 

<  l 

<  1 

20 

0 

- 

- 

- 

- 

- 

- 

_ 

■ 

Table  3,  Number  of  Drops  (Liquid  or  Equivalent  Liquid)  per  n 
or  Exceeded  With  0.  5  Percent  Probability,  Worst  Month,  Wet 


Rainfall 
Intensity 
(mm /min) 
0.60 


Median 

Diameter 

(mm) 


i  Diameter  (mm) 

0.  5-1,4 

1,5-2.  4 

2.  5-3,  4  3,  6-4.4 

4. 5-5, 4 

5, 

l-G.  4 

2020 

342 

46  6 

1 

< 

2061 

372 

48  6 

l 

< 

3062 

399 

62  7 

1 

< 

3298 

429 

66  7 

1 

< 

3634 

473 

62  8 

1 

< 

2526 

329 

43  6 

1 

< 

1650 

216 

28  4 

c  J 

< 

1002 

116 

13  2 

<  1 

< 

624 

mm 

8  1 

<  1 

< 

260 

3  <  1 

<  1 

< 

41 

<1  <1 

<  1 

< 

- 

m 

' 
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The  dashed  curves  in  Figure  1  are  based  upon  Panama  data  only}  therefore, 
probability  percentages  on  each  curve  represent  only  that  area  of  the  moist 
tropics.  Intensities  for'  these  probabilities  are  lower  than  those  for  world-wide 
rainier  tropical  areas  derived  for  MIL-STD-210B  by  Salmela  et  al  (1971)  and 
Used  by  Sissenwine  (1972).  Comparison  of  the  MIL-STD-210B  1.0,  0.5,  0. 1 
percent  and  world -record  intensities  for  the  rainiest  month  to  Panama  intensities 


(mm/hr)  specified  in  Figure  3 

L  follows: 

Probability  (%) 

1.0 

0.  5 

0.1  0.001 

World  Record 

Panama 

8.4 

37 

84  132 

— - 

MIL-STD-2  10B 

30 

48 

188 

1872 

A  0.  001  percent  world-wide  worst  month  value  was  not  detex  mined  for  MIL-STD- 
2 10B. 

Drop  numbers  in  Tables  1,  2  and  3  are  provided  for  diameters  up  to  6.4  mm, 
but  drop  diameters  greater  than  5.  5  mm  are  unstable  according  to  Fletcher  (1962) 
Curves  in  Figure  1,  however,  fail  to  indicate  a  sharp  drop  in  numbers  of  the  very 
large  sizes.  This  discrepancy  is  probably  a  result  of  an  inadequate  sample  size. 


3.  LIQUID  WATER  CONTENT 

Several  equations  have  been  developed,  mostly  by  radar  meteorologists,  for 
relating  liquid  water  content  to  precipitation  intensity.  Four  of  these  are 
presented  by  Sissenwine  (1972).  The  relationship  derived  from  Eq.  (4),  Jones 
(1956)— which  Sissenwine  found  best  supported  in  the  scientific  literature  for 
these  very  severe  thunderstorms— is  Eq.  (6).  Approximations  involving  the  cut¬ 
off  of  maximum  and  minimum  drop  sizes  were  made  to  obtain  the  drop-size  dis¬ 
tributions  from  Eqs.  (4)  and  (5)  for  Tables  1,  2  and  3.  The  importance  of  these 
approximations  can  be  estimated  from  a  comparison  of  precipitation  mass  at  the 
surface  in  the  three  storms  of  Tables  1,  2  and  3,  obtained  by  adding  the  volume 
of  the  raindrops,  to  the  values  from  Eq.  (4): 


Liquid  Water  Content 


0.  5% 

0.  1% 

1-mln 

World  Record 

(g/ni  ) 

(g/m  ) 

<g/tn3) 

Equation  (4) 

2.  2 

8.  3 

77.  2 

Drop  mass  total 

1.  8 

8.  7 

169.  0 

The  values  for  the  water  mass  from  the  drop  total  were  obtained  by  multi¬ 
plying  the  number  of  drops  in  a  0,  1-mm  diameter  range  by  the  volume  of  a  drop 
with  a  diameter  in  the  center  of  that  range,  and  summing  the  total  within  the  cutoff 
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dtamotors;  that  is,  0.  5  to  0.4  mm.  Tim  water  mass  obtained  using  this  procedure 
should  be  less  than  the  amount  obtained  using  Eq,  (6),  alnoe  only  a  certain  range 
of  drop  sizes  was  used  in  this  computation  whereas  the  entire  spectrum  of  drop 
sizes  is  inherent  in  using  Eq.  <6 >.  Since  there  Ib  a  logarithmic  increase  in  the 
number  of  drops  with  decreasing  diameter,  howovor,  there  is  a  bias  in  our  method 
of  computing  the  water  mass.  The  reason  for  the  bias  is  our  use  of  the  diameter 
in  tho  conter  of  the  range  to  compute  the  drop  volume  for  all  the  drops  in  that 
x’ango.  The  result  is  an  overestimate  of  the  drop  volume  with  Increasing  rainfall 
rate.  This  is  apparent  in  the  above  comparison  especially  for  the  1-min  world 
record.  The  discrepancy  at  the  more  realistic  extreme  rainfall  intensities  is 
small.  With  tills  in  mind,  it  can  be  seen  that  the  difference  between  water  mass 
computed  using  Eq.  (6)  to  that  using  the  drop  total  is  not  of  practical  importance. 


4.  CONCLUSIONS 

The  distribution  of  drop  sizes  in  Tables  2  and  3  gives  reasonable  numbers  of 
drops  for  the  spectrum  of  diameters  in  the  "0, 1  and  0.  5  percent  Worst  Month 
Tropics"  synthetic  storms  proposed  for  MIL-STD-2  10B  (Sissenwine,  1972). 

The  0.  b  percent  intensity  storm  was  suggested  as  mandatory  for  operation  of 
combat  equipment.  The  values  in  Table  1  give  very  improbable  all-time  extremes 
that  can  be  considered  a  "ball  park"  goal  when  life  would  be  endangered  by  an 
encounter  with  exceptional  extremes.  The  values  in  Table  2  are  for  as  Intense  a 
storm  as  could  be  assigned  a  reasonably  reliable  probability.  They  certainly 
should  be  attained  as  design  criteria  for  equipment  designed  for  world-wide  usage 
in  which  failure  would  endanger  life. 
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